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EFFECT OF VIEWING ANGLE 



PARTICLE RADIUS, 



effect of spectral range 



PARTICLE RADIUS, r(fi.m) 

FIG. 3 


EFFECT OF SPECTRAL RESOLUTION 
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PARTICLE RADIUS 



SPECIFICATIOi: 


TO ALL L’iiOM IT MAY CO’.’CERIi : 


BE IT 

Morocco, residin 
California 91106 


Californ! 


ha3 


KL'OW'i THAT ALAIN L. FYM/M , a citizen of 
g at 1080 South El Molino Avenue, Pasadena 
, in the County of Los Angeles , State of 
invented a new and useful 


FRLqUELCY-SCAL'NIMS 
SIZE SFEC7R0ME 


PARTICLE 

TER 



of which the following is a specification: 


REPi: 

ORIC.-. 


A - 


ME DISCLCiURZ 


- — ii/M 
■ ih’ POOR 




A particle size spectrometer having a fixed fiilc of 
view within the forward light scattering cone at an angle 9g 
between approximately IQQ and 200 minutes of arc (preferablv at 
150 minutes), a spectral range extending approximately froi:> 0.2 
to 4.0 inverse m.icrcmeters and a specti'ai resolution 

between about 0.1 and 0.2pm~^ (preferably toward the lower end 
of this range of spectral resolution) , is er.ployed to determine 
the distrioution cf particle sizes, independently of the chemical 
composition of tlie particles, from measurem.ents of incident light, 
Iq, at each fresuency, c ( = 1/).), and scattered light, 1(c), 
according to the equation: 


"^n ( r ) = - 


J (v) Y (v)v — lu_(oIg)da ,y=S,r/ 7 , 
- da f 


2TTr ^ o 

where i = 2a3ir.0, ? being the fixed viewing angle 0s at v;hich 
scattered ligh': is measured, r is particle size, a j.s the 
reciorocal of wavelcnrth, J is a Lessel function cf first 



|L|^ 


- 1 - 


kind and order unity, is a Bessel function of second 
kind and order unity. The cuantity, 1^, is the ratio of 
scattered light to incident light at each frequency interval. 
The apparatus is a passive rer.ote senscr that can be used in 
laboratories, field stations, flying aircrafts and airships, 
and on board an orbiting satellite. 


I 

I 


ORIGi:. OF THL IliVLMIOM 


The invention described herein '..’as r.ade in the 
perfornance of work under a NASA contract and is subject to 
the provisions of Section 3C5 of the National Aeronautics and 
Space Act of 195£, Public Law 85-558 ( 72 Stat. 435 ; USC 2457). 


BACKGROUND OF THE IKVENTIO:: 


This invention relates to a riethod and apparatus for 
determining the distribution of particle sizes of a scattering 
medium of unknown chemical composition. The medium m.ay consist 
of a composite mixture of arbitrary num.bers of chemically 
different particles eribedded in various gases. 


There is a continuing need to deter.mine the size 
distribution of any arbitrary number of particle species 
v.’hich may present a distribution in their refractive indices , 
and which m.ay be embedded in a scattering m.edium of unknown 
chemical com.position , for such industrial purposes as environ- 
mental pollution monitoring and control, chemical analyses, 
aeronautical engineering and rocket engineering as for 

such other purposes as meteorological resaarch, geophysical 
research, anc bior.edical analyses. The composite mixture of 


particles and gases may be contained in an experimental chamber, 
suspended freely in the air, or suspended in fluids such as 
medical plasmas, or in any other similar situation where parti- 
culate sizes must be determined without disturbing the sample. 


5 


10 


Conventional direct sampling methods wh^e particles 
^ are thermally precipitated, impacted or otherwise mechanically 

* collected, as well as in situ imaging methods, present several 

^ ■■■ 

^problems. Some of these problems are: disturbing of sample, 

low sampling rate, discontinuous sampling, collection efficiency 
problems, long sampling time and painstaking analyses. 

Additionally, they cannot be used as remote sensors for real- 

.. . . ' 

time operation. Oi-’ 77//.' 

^ - K POOR 


To avoid these problems of direct sampling, and at 
the same time provide a remote sensing capability, methods 
15 have been developed which exploit the forward light scattering 

properties of particles at frequencies within the range 0.193 to 
4.070 inverse micrometers (um”^) or within any larger range 
extending down to 0.02pm“^ and up to 9.27pm Incident light, 
I^, possessing a spectrum of frequencies in the ranges above 
20 specified, is characterized only by its wavenumber, Oj: , the 

reciprocal of the wavelength, X^. Its direction of propagation, 
0£, is the direction of reference, 0^=0 . The subscript, i, 
here stands for "incidence." Particles are characterized by 
their nature and micro-structure, respectively represented by 
25 a refractive index, m,.a complex number, and a .^ize distribution, 
j^n (r), representing a partial concentration per unit volume and 

i. 

per unit increment of the radius r. Only n (r) is independent 


- 3 - 


of th€ wmv«numbar. Th« various par'^lcles, assumsd to ba 
sphsrical in shapa, may also prasant a distribution in m itself. 
Gasas will ba similarly characterized by their number density, 
i.a. number of atoms or molecules yer unit volume, refractive 
index and geometrical shape. Scattered light, on the other 
hand, is characterized by its wavenumber, Cg, and (fraction, 

0g» of propagation. Here the subscript, s, stands or 
If “scat taring. " l.^e direction, 8g , is th’ s the scattering angle 
^with reference to the incident light direction 8^. 

For a given incident light, i.e., for a situation 

where and 8^ are given, «ind for a given assembly of 

scatterers, i.e., for an assembly of particles of given 

refractive indices, m(a ), aind size distribution, n (r), Mie's 

theory of scattering oy a conducting sphere enables us to 

determine the scattered light for any arbitrary Og and 8^. 

Conversely if the incident light is known for any and 8^ 

to be prescribed, and the scattered light is measured for various 

Og and 8g to be also prescribed, we should be able in principle 

to infer the properties of the particles, that is their size 

distribution, and refractive index values at the various u_. 

In other words, Mie's theory enables us to compute the output, 

I , in all its detailed variations if we know the input, I , and 
s’ o 

the properties of the particles mCCg) andn(r). However, 
what is of interest is a determination of the scatterers 
from a knowledge of the input, , and the outout , I_. But 

from an inversion of Mie's theory it would be possible to 
determine only the combined relation between the refractive 
indices, m(Og), and size distribution, n(r), and noj^ the 
size distribution separatel''. 


It has b««n discovarad that the Fraunhofer theory of 


10 


15 


angular diffraction of light of fixed single frequency at an 

aperture in a plane screen can be extended to a range of multiple 

frequencies. Under these conditions, the light scattered by a 

particle is essentially independent of the refractive index, 

mCOg), and depends only on its size as though the particle 

Ipere an aperture in a screen of the same radius. It has also 

^een discovered that, working within the restricted domain 

of applicability of the theory thus extended (9 between 100 and 

200 minutes of arc, approximately, preferably 9g * 150 minutes; 

emd 0.2pm“^ 1 Og 1 4.0um“^, approximately), it is possible to 

determine size distribution, n (r), independently of refractive 

indices m(a ), i.e., to determine size distribution of particles 
s 

of radius larger than approximately one micrometer (um) without 
knowing anything of their refractive indices, from known incident 
light, Iq, and measured scattered light, Ig. 


20 


If the assembly of ccatterers under study consists of 
particles of j different species and k different gases, the 
identical determination of n (r) can still be made provided only 
j, k, and the so-called depolarization factor of each gas (this 
is defined as the ratio of scattered intensities in directions 
parallel and perpendicular to the plane of scattering for an 
incident beam of natural light), Pj^, are known. 


SUmARY OF THE INVENTION 




i'iw 

Poor 


25 In accordance with the present invention, a natural 

or artificial light source is passed through a volufde of 

scattering particles the size distribution of which is to be 
determined. The intensity of light scattered by the particles. 


at a fixed angle between approximately 100 and 200 minutes of 
arc from the direction of incident light, preferably 150 minutes 
of arc, is recorded as a function of frequency, 1(a). Records 
are made across the approximate frequency interval 0.20um”^ 
to 4.0um“'^, at frequency steps smaller than about 0.2um"^. 
Simultaneously, the intensity of the incident light at each 
^frequency, Iq, is recorded. The ratio of scattered light, 
^l(a), to incident light, 1^, is then obtained. In order to 
determine the size distribution, n (r), fixim the ratio 1^^, the 
following eouation has been derived and investigated with regard 
to its correctness, the uniqueness of the solution obtained from 
it, and the stability of this solution with regard to inherent 


experimental and analytical noises '• 


i'U> ilv 


n(r) = - 


Jl(y)Yi(y)y-~ (ol^)do, y=£ro. 


where is a- Bessel function of first kind and order unity, 
is a • Bessel function of second kind and order unity, 

^ - 2TTsin9, 0 being the fixed viewing angle at which scattered 
light is measured, r is particle radius, and a is the reciprocal 
of wavelength. In the case of j particle species embedded in k 
gas species , the former equation remains valid provided only 
is changed to the following expression: 


1 r4rnc;0 T - 1 V' ^ “ Pk COS28 

tL " 2 * - ,; 


where Pj^ is a so-called depolarization factor of the k^“ gas 
defined as the ratio of scattered intensities in directions 
parallel and perpendicular to the plane of scattering for an 
incident beam of -natural light. * 


1 ^ 
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The novel features of the invention are set forth 
with particularity in the appended claims. The invention 
will best be understood from the following description when 
read in conjunction with the accompanying drawings. 

% 

% 

5 . BRIEF DESCRIPTION OF THE DRAWINGS 

■ I — »ii I 

The FIG. 1 illustrates schematically an example of 
apparatus for the method of the present invention, where the 
method steps operating on measu'.'ed quantities are represented 
by functional blocks. 

10 FIG. 2 is a graph showing the effect of viewing 

angle, 0 , on determining particle size distribution from 
light scattered by the particles according to the present 
invention . 

FIG. 3 is a graph showing the effect of spectral 
15 range to Q^iax* scanning the frequency of light 

incident on particles for determining particle size distri- 
bution. 

FIG. 4 is a graph showing the effect of spectral 
resolution, Aa , in the process of determining particle size 
20 distribution according to the present invention. 


- 7 - 


DESCRIPTION OF THE PREFERRED EMBODIMENTS 


Referring to FIG. 1, a natural or artifi?J.al 
light source 10 is directed to a cell 11 filled with a medium 
of unknown chemical composition through a chopper 12 which 
cyclically deflects the incident light, I^, to a reference 
^ multiple-filter detector 13 for the purpose of measuring the 
incident light intensity as the light source is scanned in 
frequency. The filter of the detector is preferably changed 
as the light source is scanned in order that the reference 
measurement I^ be accurate for each frequency of incident 
light. An object multiple-filter detector 14 is likewise 
controlled to detect scattered light at the different 
prescribed frequencies as the source is scanned. That is 
represented by dotted lines 15 and 16 for synchronized control 
of the detector and light source. The reference and object 
multiple-filter detectors can be detector arrays in which each 
detector element is slaved to a specific frequency. 


Ay- 


' The iql^^^Adetector 14 is designed to have a fixed 


' ’■ 1 


field of view at a fixed angle 0 of less than 6® within the 

forward scattering cone of the cell 11. For best results, the 

scattering angle 9 is made about 150 minutes of arc, but it 

s 

may be any angle between 100 and 200 minutes of arc, approxi- 
mately. It is also designed to have a spectral resolution, 
i.e., a frequency step, Aa, as coarse as 0.2um"’^ or as fine as 

O.lym”^, and to operate within the approximate spectral frequency 

-1 -1 
limits of a = 0.2pm and “ 4.0pm 

mxn ina.x 




1 


Although a cell is indicated, the scattering particles 


forming a homogeneous cloud or medium between the light source 
and the detector may be contained in a laboratory chamber, 
suspended freely in the air, contained in fluids such as plasmas 
5 ^d blood, emitted from any type of industrial or engine 

nozzle, or in any other similar situation where particulate 
^ sizes must be determined without disturbing the sample in any 
manner. Also, although a variable light source is indicated, 

' radiance from particles that exhibit significant thermal emission 
10 e.g., radiance from atmospheric particles which emerges from 

the atmosphere, can also be used. On the other hand, a light 
source may be employed that emits all of the frequencies to be 
scanned. In the latter case, the control for frequency scanning 


may be effected at only the detectors 13 and 14, while in the 


former case it is effected only at detector 14. 
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If the foregoing restrictions on the fixed angle, 


6g, spectral frequency range, <^niin ^max» spectral resolu- 
tion, Ao, of the reference and object detectors are observed, 
the measurements of scattered light 1(a) can provide the complete 
20 spectrum of particulate sizes, i.e., particle size distribution, 
n(r), in radii of equivalent spheres, from approximately lym 
to any arbitrarily larger value, including the accurate location 
of the mode radius (corresponding to the largest concentration 
of particles). The technique will now be described. Consider 
2 5 an inciden .* beam of radiation of wavenumber, a = 1/X, (X=wave- 
length), and intensity, 1^, illuminating a single particle of 
radius, r. The singly-scattered light intensity, I(o), at 
angle, 0, in units of 1^, is expressed by the formula: 
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r Ji Ckr »ln 6) 1 2^ 
sin 8 J 

where k = 2ita and is a Bessel function of first kind and 
order iinity. For a distribution of sizes, n(r), the generalize 
f tion of this expression will provide the total scattered light 
5 intensity : 
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n(r) r J^(x)dr, x=kr sin 0. 
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It may be noted that equations (1) and (2) depend only on r. 

In order to determine n(r) from measurements of I^, it is 
necessary to invert this equation. The analytical inverse of 
10 equation (2) is obtained in the form: 

3 ^ eo 

n(r) = - - J: -..- J y=Jlra, ( 


where Jl = 2tt sin 0 and is a . Bessel function of second kind 
and order unity. 
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If the scattering medium being illuminated consists 
of j different species of particles embedded in k different 
species of gases, equations (2) and (3) remain valid provided 
only we replace in these equations the quantity I^ by the 
following quantity: 
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G„ = ± I 4 cos 0 I^ - ■J 
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where p, is the depolarization factor of the k gas 
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It must be stressed that Equation (3) utilizing either 
I- is a closed form expression for n(r) which makes no 

assumption on any amalytical form by which this function nay 
be modeled. It is exactly true for the function I defined by 
Equation (2) or the function defined by Equation (4). It 

may be noted, howeyer, that a will not yary oyer an infinite 
range as indicated by Equation (3). ^ OOR 

For example, if it is wished to use this instrument 
as a remote sensor on board a planetary orbiter, o has a 
lower bound which will yary according to the brightness temper- 
ature of the planetary atmosphere of interest. This yalue of 
o is that at which the planetary radiation dominates the yisible 
solar radiation reflected by the planet. This is because the 
thermal radiation from the planet surface and atmosphere system, 
which is emitted due to the temperature of this surface and 
atmosphere, is not the absolute zero. Thus, on Earth, this cut- 
off yalue is approximately 0.37pm Howeyer, although the 
integration indicated in Equation (3) cannot be carried out 
completely, and it is conceiyable that in some extreme cases 
the distribution n(r) may either not be reproduced exactly or 
not at all, it has been discoyered that an object detector 
designed to operate within the limits set forth hereinbefore 
will permit the integration to be carried out to proyide the 
detailed spectrum of sizes (particle size distribution) in 
equiyalent spheres of the scattering particles. These limits 
proyide a unique method which has been submitted to detailed 
scrutiny as to stability with regard to both inherent experi- 
mental and analytical noises. 


74/348 


The first computational step of forming the ratio 
indicated by Equation (2) is represented by a functional block 
17, and the second computational step indicated by Equation 
(3) is represented by a functional block 18. The ratio can 
5 be formed using analog circuits, or upon converting the detector 
signals to digital form, by digital techniques. In either case, 
the result of the first computational step is preferably obtained 
in digital form for storage and subsequent retrieval for solu- 
"i ..on of Equation (3). The values of this ratio at the 
10 several prescribed frequencies are then inserted in Equation 

(3) which is subsequently solved for various values of r larger 
than approximately lum. For example, for r = r^^, the equation 
yields n^ = n(r=r^); for r = r 2 » the solution yields n 2 = n(r=r 2 ) 
and so forth. The set of values n^^, n^j ... is the desired 
15 particle distribution n(r). A general purpose computer 

can be used for solving Equation (3). The computer can be 
used on a time-sharing basis or a special-purpose computer 
could be slaved to the apparatus. The output from such a 
computer is converted to an electrical signal which is then 
20 displayed visually on a cathode ray tube (CRT) display 19. 

Thus, a real-time visualization of the size distribution is 
provided, such as for the purpose of monitoring the changes in 
particle sizes that may take place as is the case, for example, 
with environmental particulate pollution. 

25 The effect of viewing angle, 6^ , can be appreciated 

from the graph of FIG. 2 which shows a true plot of particle 
size distribution by a bell-shaped curved labeled "n(true)” 
for comparison with particle size distribution curves determined 
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at viewing angles of 50 minutes, 150 minutes and 350 minutes. 

The distribution curve labeled "0g = 15 0 min’* is found to be 
very close to the true distribution curve. Consequently, it 
is concluded that a viewing angle of about 150 min is optimum. 

At viewing angles between approximately 100 min and 200 min, 
the results are still very good, although small negative tails 
develop close to these two values below 1.5um radius; but 
beyond that range there is not only further degradation of the 
results from the true curve, but totally incorrect results over 
much of the range of particle sizes. For example, at a viewing 
angle as small as 50 r i.n, the results indicate a negative 
distribution between 14 and 17um radius, and although it seems 
to follow the true distribution curve as particle radius 
decreases, it becomes clear at below Bum radius that the distri- 
bution curve is totally in error (meaningless) since size 
distribution increases, apparently without limit, as the particle 
radius decreases. Deviation from the optimum viewing angle in 
the other direction will not produce such absurd results, as 
shown by the distribution curve labeled "6 = 350 min". But 

still the results have so departed from the true distribution 
as to not be useful except perhaps in locating the peak 
position. ' Note the negative distribution above Hum and below 
2.7um radius. The true distribution cannot, of course, go 
below zero. This demonstrates the limited range of viewing 


angles that can be effectively used in this technique. 
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The effect of spectral range can be appreciated from 
the graph of FIG. 3. In the range of about 0.2 to 4.0uni 
the result is very nearly identical to the true distribution 
curve as shown by a curve labeled "0.198 ^ a ^ 4.070um 
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This result remains substantially the same if the starting value 
of 0.198 were to be lowered down to 0.02iim”^. At the lower 
part of the range, namely at 0.198 to 2.860wm the form of the 
distribution curve remains about the same but becomes negative 
at leas than 2ym radius. Increasing the upper end of the range 
to approximately O.OOpm ^ produces about the same results, 
but the distribution curve becomes more negative at less than 
2pm radius. Working in a subinterval of the lower part of the 
optimum range, namely in the range from 0,44 to 2.86pm the 
result is about the same as working at just the lower part of 
that optimum range, except that the distribution is found to 
be negative below 3pm radius and above 12pm radius, with 
decrease of approximately 18% in the peak. Operating at only 
a narrow central part of the optimum spectral range, namely 
between about 0.3 to 1.2pm~^ produces the absurd result 
indicated by the curve labeled ”0.286 - a - 1.190pm 
Operating below the optimum range, such as in a range from 0.03 
to 0.2pm"'^ also produces an absurd result, as shown by the 
curve labeled ”0.031 - a ~ 0.235pm Lastly, extending 

the optimum range, for example from 1.0 to 10.0pm ^ also produces 
the absurd respll^ shown by the curve labeled ”1.0 - a - 10.0pm 


The effect of spectral resolution is illustrated in 
FIG. 4. Again the curve labeled n(true) is the theoretical 
or true particle distribution. The optimum is shown to be a 
spectral resolution of about 0.06pm as shown by the curve 
labeled "La ~ 0,061pm However, spectral resolution is less 

critical than spectral range or viewing angle. A smaller 
resolution would have no deleterious effect on the result. 


A l«rg«r r«8olution can b« tolarated to about thra* times 
the optimum. There, at the resolution of 0.18um“^, the 
results begin to degrade, and at about O.SOjjm the 
results become absurd as shown by the curve labeled 
5 "La - 0.774ym”^." 

The necessary spectral range with the optimum spectral 
resolution may be achieved by scanning frequency at the 
detector using a detector-array in which each detector is 
slaved to a specific frequency, as noted hereinbefore. Alter- 
10 natively, one or more voltage variable photodetectors are used, 
and frequency scanning is achieved by varying the bias voltage 
to the detectors. For most laboratory applications, there will 
be a specific light source, either variable in frequency or 
containing all frequencies of interest. In the latter case, 

15 scanning the frequencies in the desired range could be accomplishe 
by the use of filters. Such an arrangement would also be used 
outside the laboratory from a ground-based station such as for 
determining the size distribution of particles in the Earth’s 
atmosphere, where light from the sun is filtered for scanning 
20 the frequency of light passed through a volume of scattering 

particles. In the case of an orbiting satellite, on the other han 
it is the thermal emission from the atmosphere-surface system 
that would be filtered. 

Although particular embodiments of the invention 
25 have been described and illustrated herein, it is recognized 
that modifications and equivalents may readily occur to those 
skilled in the art and consequently it is intended that the 
claims be interpreted to cover such modifications and equiva- 
lents. 


- 15 - 


